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ABSTRACT 

Context. Since 1995, more than 1500 exoplanets have been discovered around a large diversity of host stars (from M- to A-type stars). 
Tidal dissipation in stellar convective envelopes is a key actor that shapes the orbital architecture of short-period systems. 

Aims. Our objective is to understand and evaluate how tidal dissipation in the convective envelope of low-mass stars (from M to F 
types) depends on their mass, evolutionary stage and rotation. 

Methods. Using a simplified two-layer assumption, we compute analytically the frequency-averaged tidal dissipation in their con¬ 
vective envelope. This dissipation is due to the conversion into heat of the kinetic energy of tidal non wave-like/equilibrium flow and 
inertial waves because of the viscous friction applied by turbulent convection. Using grids of stellar models allows us to study the 
variation of the dissipation as a function of stellar mass and age on the Pre-Main-Sequence and on the Main-Sequence for stars with 
masses spanning from 0.4 to 1.4Mq. 

Results. During their Pre-Main-Sequence, all low-mass stars have an increase of the frequency-averaged tidal dissipation for a fixed 
angular velocity in their convective envelope until they reach a critical aspect and mass ratios (respectively a = Rc/Rs and j8 = MdM^ 
where R^,M^,Rc and are the star’s radius and mass and its radiative core’s top radius and mass). Next, the dissipation evolves 
on the Main Sequence to an asymptotic value that becomes maximum for 0.6Mq K-type stars and that decreases by several orders 
of magnitude with increasing stellar mass. Finally, the rotational evolution of low-mass stars strengthens the importance of tidal 
dissipation during the Pre-Main-Sequence for star-planet and multiple star systems. 

Conclusions. As shown by observations, tidal dissipation in stars varies over several orders of magnitude as a function of stellar 
mass, age and rotation. We demonstrate that: i) it reaches a maximum value on the Pre-Main-Sequence for all stellar masses; ii) on 
the Main-Sequence and at fixed angular velocity, it becomes maximum for O. 6 M 0 K-type stars and it decreases with increasing mass. 
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1. Introduction and context 

Since 1995, more than 1500 exoplanets have been discovered 
around a large diversity of host stars (e.g. |Perryman||201 l| l. 
As of today, their mass range spreads from M red dwarfs to 
intermediate-mass A-type stars. In this context, tidal dissipation 
in the host star has a strong impact on the orbital configuration 
of short-period systems. Indeed, the dissipation in rotating tur¬ 
bulent convective envelopes of low-mass stars is believed to play 
a key-role for tidal migration, orbits’ circularization and spins’ 
alignment (e.g. [Albredit et al.|[2012[ |Lai||2012{ |Ogilvie|[2014[ 
and references therein for hot-Jupiter systems). In these regions, 
this is the turbulent friction acting on tidal flows, which dissi¬ 
pates their kinetic energy into heat, that drives the tidal evolution 
of star-planet systems. In stellar convection zones, tidal flows 
are constituted by a large-scale non-wave like/equilibrium flow 
driven by the hydrostatic adjustment of the stellar structure be¬ 
cause of the presence of the planetary/stellar companion ( Zahn| 
|1966[ [Remus et al. 12012} and the dynamical tide constituted by 
inertial waves, their restoring force being the Coriolis accelera¬ 
tion (e.g. Ogilvie & Lin|2007 1. In this framework, both the struc¬ 
ture and rotation of stars strongly varies along their evolution 
(e.g. Siess et al.|2000| Gallet & Bouvier|2013 20151. Moreover, 
as reported by Ogilvie ( 2014| l, observations of star-planet and bi¬ 
nary star systems show that tidal dissipation varies over several 


orders of magnitude. Therefore, one of the key questions that 
must be addressed is; how does tidal dissipation in the convec¬ 
tive envelope of low-mass stars varies as a function of stellar 
mass, evolutionary stage and rotation? 

In this work, we study the variations of the frequency- 
averaged tidal dissipation in stellar convective envelopes as a 
function of the mass, the age and the rotation of stars. In sec. 
2, we introduce the assumptions and the formalism that allows 
us to analytically evaluate this quantity as a function of the struc¬ 
ture and rotation of stars ( Ogilvie|2013| l. In sec. 3, we compute 
it as a function of stellar mass and evolutionary stage at fixed an¬ 
gular velocity using grids of stellar models for stars from 0.4 to 
1.4Mo. Next, we discuss the impact of the rotational evolution of 
stars. In sec. 4, we present our conclusions and the perspectives 
of this work. 


2. Tidal dissipation modelling 

To analytically evaluate the frequency-averaged tidal dissipation 
in the convective envelope of low-mass stars, we here choose to 
adopt a simplified two-layer model as in Ogilvie ( |2013| l (see also 
jPenev et al.|2014| and fig.[2l- In this modeling, both the radiative 
core and the convective envelope are assumed to be homoge¬ 
neous with respective constant densities pc and pe for the sake of 
simplicity. This model features a star A of mass Ms and mean ra- 
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Fig. 2. Top-left: track of evolution of stars from 0.4 to 1.4Mq in the Hertzsprung-Russell diagram that gives the luminosity (L) as a function of 
effective temperature (Teff) (red, orange, green, dark blue, cyan lines correspond to M, K, G, F and A-type stars respectively). Top-right: Evolution 
of the stellar radius of stars from 0.4 to 1.4Mq as a function of time. Bottom-left: Evolution of the radius aspect ratio a = Rc/Rs of stars from 
0.4 to 1.4Mq as a function of time. Bottom-right: Evolution of the mass aspect ratio p = M^jM^ of stars from 0.4 to 1.4Mo as a function of time. 



Fig. 1. Two-layer low-mass star A of mass and mean radius R^ and 
point-mass tidal perturher B of mass m orbiting with a mean motion n. 
The radiative core of radius Rc, mass Me and density pc is surrounded 
by the convective envelope of density Pe. 


dius Rs hosting a point-mass tidal perturher B of mass m orbiting 
with a mean motion n. The convective envelope of A is assumed 
to be in moderate solid-body rotation with an angular velocity 


Q so that 




« 1 


where Qp 


the critical angular velocity and Q is the gravitational constant. 
It surrounds the radiative core of radius and mass M^. 

Tidal dissipation in the external convection zone of A origi¬ 
nates from the excitation by B of inertial waves, which have the 
Coriolis acceleration as restoring force. They are damped by the 
turbulent friction, which is modeled using a turbulent viscosity 
(see|Ogilvie & Lesur|2012[ and references therein). Its analyti¬ 


cal evaluation in our two-layer model was conducted by Ogilvie 
( |2013| l who assumed an homogeneous and incompressible con¬ 
vective envelop^ surrounding an homogeneous stable fluid core 
where no dissipation occurs. The solutions of the system of dy¬ 
namical equations for the envelope written in the co-rotating 
frame are separated into a non-wavelike part (with subscripts 
nw), which corresponds to the immediate hydrostatic adjustment 
to the external tidal potential (t/), and a wavelike part (with sub- 


* In this regime, the centrifugal acceleration, which scales as is 
neglected. 

^ This assumption corresponds to inertial waves with shorter wave¬ 
length than the characteristic length of the variation of density. 
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script „) driven by the action of the Coriolis acceleration on the 
non-wavelike part (see Appendix [A|and [Ogilvie ( 2013| l). The ki¬ 
netic energy of the wavelike part of the solution can be derived 
without solving the whole system of equations, thanks to an im¬ 
pulsive calculation. It is dissipated after a finite time and allows 
us to compute the frequency-averaged tidal dissipation given in 
Ogilvie| ( |2013t (eq. B3f] 
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We introduce the tidal frequency u - sn - inQ. (with s eZ) and 
the Love number k'”, associated to the {l,m) component of the 
time-dependent tidal potential U that corresponds to the spheri¬ 
cal harmonic Y'". It quantifies at the surface of the star A (r = R^) 
the ratio of the tidal perturbation of its self-gravity potential over 
the tidal potential. In the case of a dissipative fluid as in convec¬ 
tive envelopes, it is a complex quantity that depends on the tidal 
frequency with a real part that accounts for the energy stored in 
the tidal perturbati on while the imagin ary part accounts for the 
energy losses (e.g. 


Remus et 


aL||2012 1. Note that Im|^k™(w)j is 


proportional to sgn(w) and can be expressed in terms of the tidal 
quality factor Q’l'ico) or equivalently the tidal angle dfito), which 
both depend on the tidal frequency; 


Qfioj) ^ = sgn(w) \kJ‘(aj)\ * Im [k)"(w)] = sin [2 df (w)]. (3) 

^ In the calculation, the weighted frequency integral reduces to w 6 
[-2n, 2n] because higher-frequency acoustic waves are filtered out. 
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We here choose to consider the simplest case of a coplanar sys¬ 
tem for which the tidal potential (U) reduces to the component 
(/ - 2,m - 2) as well as the quadrupolar response of A. To un¬ 
ravel the impact of the variation of stellar internal structure, we 
introduce a first frequency-averaged dissipation at fixed normal¬ 
ized angular velocity 

= e-2 (D)^ = e-2 (im [kl(aj)])^ , (4) 


which depends on a and /? only, and a second one 






-1 


R. 


Rc, 


(5) 


where = ^Q/ ^0Mq/R^ - (Mq, Rq and Q® being 

the solar mass, radius and critical angular velocity), that allows 
us to take also into account the dependences on the radius (Rs) 
variations and on the mass {Ms). Note that these quantities being 
averaged in frequency, the complicated frequency-dependence 
of the dissipation in a spherical shell ( jOgilvie & Lin|[2007| l is 
filtered out. The dissipation at a given frequency could thus be 
larger or smaller than its averaged value ((2))(j) by several orders 
of magnitude. Moreover, constitutes a lower bound of the 
typical dissipation in the whole star since it is also necessary to 
take into account the damping of tidal gravito-inertial waves in 
radiation zones ([Ivanov et al.|2013|l. 



-2 

-3 

-4 

-5 

-6 

-7 

-8 


Fig. 3. Variation of (2)>“ = e ^ (im as a function of aspect and 

mass ratios {a and fi respectively) in color scales. Evolutionary tracks 
of stars from 0.4 to l.TMg in the {a,fi) plane (the colors are the same as 
in Fig|^. 


mass along the PMS. Next, for stars with Ms > 0.5Mq, both a 
and yS reach a value on the Zero-Age-Main-Sequence (ZAMS) 
that stays almost constant along the MS. For example, in the 
case of a IMq star we recover the usual solar values a ^ 0.71 
and p 0.98. Both a and p increase with Ms on the MS. Finally, 
when stars reach their SG phase, both a and p decrease because 
of the simultaneous extension of their convective envelope and 
contraction of their radiative core. 


3. Tidal dissipation along stellar evolution 

3.1. Grids of models for low-mass stars 

To study the variation of tidal dissipation in the convective en¬ 
velope of low-mass stars along their evolution in the theoretical 
framework presented above we need to compute the stellar ra¬ 
dius (Rs) and the mass and aspect ratios (respectively a and P) 
as functions of time. We use grids of stellar models from 0.4 
to 1.4Mo for a metallicity Z = 0.02 computed by jSiess et ak] 
( 2000| l using the STAREVOL code. These stellar models solve, 
at the zeroth-order, the hydrostatic and energetic spherical bal¬ 
ances using the relevant equation of state and description of stel¬ 
lar atmospheres (a complete description of the STAREVOL code 
is given in jSiess et al.|200()| l. They allow us to follow the evolu¬ 
tion of stars and of their structure along their evolutionary tracks 
in the Hetzsprung-Russell diagram (Eig. top-left panel). In 
this work, we choose to focus on the Pre-Main-Sequence, Main- 
Sequence and Sub-Giant phases of evolution (hereafter reported 
as PMS, MS and SG). In Pig. (top-right panel), the variation 
of the radius of the star (Rs) is plotted as a function of time. Por 
every stars, we identify the three well-known phases: i) during 
the PMS, stars are in contraction and their radius (Rs) decreases; 
ii) on the MS it reaches an almost constant value until the begin¬ 
ning of the SG phase; iii) on the SG phase, stars’ envelope begin 
to expand and Rs increases. As this has been seen in Sec. 2., the 
important structural properties to compute for tidal dissipation 
are also the mass and aspect ratios relative to the radiative core 
on which tidal inertial waves propagating in the convective en¬ 
velope can reflect, leading potentially to wave attractors where 
strong dissipation may occur (e.g. jOgilvie & Lin||2007[ [Good-j 
man & Lackner|2()09 1 . In Pig. we thus plot the temporal evo¬ 
lution of a and p for the different stellar masses studied here (re¬ 
spectively in bottom-left and bottom-right panels). After a first 
phase where the star is completely convective on the Hayashi 
phase, the central radiative core is growing both in radius and 


3.2. Time evolution of tidal dissipation as a function of stellar 
mass and evolutionary stage 

The phases of the structural evolution of low-mass stars being 
recalled, we now have to compute the corresponding evolution 
of the dissipation of the kinetic energy of tidal inertial waves in 
their convective envelope as a function of time. In Pig.|^ we thus 
represent: i) the intensity of the first frequency-averaged dissi¬ 
pation at fixed normalized angular velocity {{D}^; see Eq.|^ as 
a function of a and p only in color scales and ii) the position of 
the couples (a,/?) along the evolution of stars for stellar masses 
from 0.4 to 1.4Mo. On one hand, (fD)^ presents a maximum 
in an island region around (amax ~ O.STTjSmax ~ 0.501). On 
the other hand, the stars of different masses have a different 
evolution of their internal tidal dissipation in their envelope 
because of their different structural evolution. Pirst, all stars 
have an increase of (fD)^ during the beginning of their PMS 
until they reach a position in the diagram close to (amax.ySmax) 
because of the formation and the growth of their radiative core. 
Then, (2))^ has an evolution that is a direct function of stellar 
mass. Por 0.4Mo M-type stars, the radiative core has an evo¬ 
lution where it Anally disappears. Therefore, the configuration 
converges on the MS to the case of a fully convective star with 
the weak dissipation of normal inertial modes ( |Wu||2005) . Por 
stars from Q.SMq to 1.4Mo masses, the evolution is different. 
After they have reached their position closest to (amax,jSmax) in 
the diagram, they evolve to a position (aMS^y^Ms) where they 
stay during almost the whole MS because of the weak evolution 
of these quantities. Since aMS Pms are increasing functions 
of the stellar mass, (29)^, which is almost constant on the MS, 
decreases from K to A-type stars. However, the evolution of Rs 
must also been taken into account. Therefore, the evolution of 

(see Eq. as a function of time (left-panel) and of the 
effective temperature (right-panel) are given in Pig. [^taking into 
account the simultaneous variations of a, p, and Rs. We recover 
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Fig. 4. Evolution of the frequency-averaged tidal dissipation at fixed normalized angular velocity, - e ^ ^Im , as a function of time 

(left panel) and effective temperature (right panel) for stellar masses (MJ from 0.4 to 1.4Mo (the colors are the same as in Fig|^. 


the two phases of the evolution of the dissipation. On the PMS, 
it increases towards the maximum value, which grows with 
stellar mass, corresponding to the region around (amaxjjSmax) 
for all stellar masses. The time coordinate of this maximum is 
smaller if stellar mass is higher because of the corresponding 
shorter life-time of the star. Then, the dissipation decreases to 
reach rapidly its almost constant value on the MS. As already 
pointed out before and expected from observational constraints 
(e.g. [Albrecht et al. [201 ^ and previous the oretical works 


(Ogilvie & Lin 2007 iBarker & Ogilvie 2009| p 


it decreases 

with stellar mass from Q.6Mq K-type to A-type stars by several 
orders of magnitude (a; 3 between 0.6 and I.TMq) because of 
the variation of the thickness of the convective envelope, which 
becomes thinner. Finally, for higher-mass stars, we can see a 


final increase 




because of the simultaneous extension of 


their convective envelope and contraction of their radiative core 
during their SG phase leading again the star towards the region 
of maximum dissipation in the ia,/3) plane. 


3.3. The impact of the rotational evolution of stars 

It is also crucial to discuss consequences of the rotational evo¬ 
lution of low-mass stars. As we know from observational works 
and related modeling (e.g. [Gallet & Bouvier 2013 2015[ and 
references therein), their rotation follows three main phases of 
evolution: i) first, stars are trapped in co-rotation with the sur¬ 
rounding circumstellar disk; ii) next, because of the contraction 
of stars on the PMS (e.g. Fig. ^ top-left panel) their rotation 
(and thus e) increases; iii) finally, stars are braked on the MS 
because of the torque applied by pressure-driven stellar winds 
(e.g. [Reville et al.|2015| and references therein) and e decreases. 
As a detailed computation of stellar rotating models is out of 
the scope of the present work because of the complex internal 
and external mechanisms that must be taken into account, we 
will not here evaluate the variation of tidal dissipation because of 
the simultaneous evolution of stellar internal structure and rota¬ 
tion. However, from results obtained by |Gallet & Bouvier| ( |2013[ 
2015| l, we can easily infer that {fD}^ is increased during the PMS 
because both of the growths of the radiative core and of the an¬ 
gular velocity. From Fig. 5 of |Gallet & Bouvier] ( |2015) l, we can 
see that in the case of the 0.5Mq star these growths will be simul¬ 
taneous. On the MS, (fD)^ decreases because of the evolution of 
the structure of the star and of its braking by stellar winds. As 
demonstrated before, the value of is determinant. 


In the cases of F and A-type stars, a convective core is present in 
addition to the envelope. We neglect its action assuming it corresponds 
to the regime of weak dissipation in a full sphere. 
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4. Conclusions 


All low-mass stars have an increase of the frequency-averaged 
tidal dissipation for a fixed angular velocity in their convective 
envelope for tidal frequencies lying within the range [-2Q., 2Q.] 
so that inertial waves can be excited until they reach a critical as¬ 
pect and mass ratios close to (amax^ySmax) during the PMS. Next, 
it evolves on the MS to an asymptotic value that becomes maxi¬ 
mum for Q.6Mq K-type stars and that decreases by several orders 
of magnitude with increasing stellar mass. Finally, the rotational 
evolution of low-mass stars strengthens the importance of tidal 
dissipation during the PMS for star-planet and multiple star sys¬ 
tems as pointed out before by Zahn & Bouchet (19891. 

In a near future, it would be important to improve the physics 
of dissipation in stellar convection zones by taking into account 
density stratification, differential rotation, magnetic field and 
nonlinear effects and related possible instabilities (e.g. |Baruteau| 
|& Rieutord|2013[|Schmitt|2070||Favier et al.|2014| ). Moreover, it 
will be necessary to include results obtained for the dissipation 
of tidal gravito-inertial waves in stellar radiative cores ( Ivanov] 
et akpoO]! to get a complete picture (|Guillot et al.|[^ 4|l. Fi- 


nally, it will be interesting to explore advanced phases of stellar 
evolution. 
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Appendix A: Dynamical equations for tidal flows in 
convective envelopes 

The solutions of the system of dynamical equations for the enve¬ 
lope written in the co-rotating frame are separated into a non¬ 
wavelike part (with subscripts nw), which corresponds to the 
immediate hydrostatic adjustment to the external tidal potential 
(U), and a wavelike part (with subscript „) driven by the action 
of the Coriolis acceleration on the non-wavelike part: 



Sw -t- 2 e^ X Sw — “VlTw + f, 


and W ^ ^ p'^ ^ 0, 

V ■ (pe Sw) = 0, 


Pnw “ ^ ' )Pe Snw)» 

V^d)' ^47tGp' 


(A.l) 


where s is the displacement, the unit vector along the rotation 
axis, h the specific enthalpy, and d) is the self-gravitational po¬ 
tential of A. Primed variables denote an Eulerian perturbation in 
relation to the unperturbed state with unprimed variables. Note 
that U being of the first order in tidal amplitude, it appears in 
perturbation equations. Finally, W = VTnw + = h' + (t>' + U 

while f = -2 Cj X Snw is the acceleration driving the wavelike 
part of the solution. 


Article number, page 5 of|^ 






